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We present an analytical and comprehensive description of the quantum dynamics of a microwave 
resonantly driven superconducting phase qubit coupled to a microscopic two-level system (TLS), 
covering a wide range of the external microwave field strength. Our model predicts several in- 
teresting phenomena in such an ac driven four-level bipartite system including anomalous Rabi 
oscillations, high-contrast beatings of Rabi oscillations, and extraordinary two-photon transitions. 
Our experimental results in a coupled qubit-TLS system agree quantitatively very well with the 
predictions of the theoretical model. 
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Microscopic two-level systems (TLSs) originating from 
variety of disorders are ubiquitous in solid state devices. 
In particular, the TLS plays a double-faced character 
in the implementation of quantum information proces- 
sor based on the superconducting qubits [THU]- On one 
side, the TLSs are harmful for the coherent control of 
qubit states |4l410] . Although significant progress has 
been made to reduce the TLSs in the superconducting 
qubits [TTJ [TJ] it seems however difficult to completely 
eliminate them in the foreseeable future. On the other 
side, the TLSs can play a useful role because of their rel- 
atively long coherence time p~3l417J . It has been demon- 
strated that the TLSs in the tunnel barrier of a Joseph- 
son junction can function as naturally formed qubits |13j . 
quantum memory cells { 14] , and quantum beam split- 
ters [TB]. In either case a thorough understanding of 
the quantum dynamics of the coupled qubit-TLS system 
is critical. What's more, the quantum dynamics of a 
resonantly driven coupled quantum bipartite system is 
important not only to fundamental physics [T5J Q]5] but 
also to applications such as building scalable quantum 
information processors based on qubits [IH3 HU] . How- 
ever, previous theoretical studies mainly focused on a 
narrow range of the ac driving strengths and only nu- 
merical simulations were performed, leaving the physics 
picture unclear. Consequently, the experimental investi- 
gations were incomplete although interesting two-photon 
Rabi oscillations and anomalous Rabi oscillations have 
been observed [H[HJ[55]. Here, we present an analytical 
and comprehensive description of the quantum dynamics 
of a resonantly driven qubit-TLS system that covers a 
wide range of ac field and qubit-TLS coupling strengths. 



More importantly, we have experimentally investigated a 
superconducting phase qubit coupled to the TLSs with 
all critical system parameters calibrated independently 
and the data strongly support the results of our theoret- 
ical analysis. Therefore, the results reported here, both 
theoretical and experimental, provide many insights into 
the dynamics of resonantly driven coupled four-level bi- 
partite quantum systems. 

We start by modeling the qubit-TLS system with 
the Hamiltonian [T5J Q71 H{t) = H q (t) + H T + 

H q ^ T . The qubit's Hamiltonian is H q {t) = —^uj q al + 
Ml m cos u)ta% , where H is Planck's constant divided by 
2ir, hu! q is the energy level spacing of the qubit, f2 m /27r 
is the microwave induced Rabi frequency, and ui is the mi- 
crowave frequency. The Hamiltonian of the TLS can be 
written as Ht = — §wtoT, where JvjJt is the energy level 
spacing of the TLS. The interaction Hamiltonian then is 
H q -T — HQ, c a% ® crj, where fl c is the coupling strength 
between the qubit and the TLS, o% yx {cr^ y z ) are the 
Pauli operators acting on the states of the qubit (TLS). 
Such o~ x coupling is usually found in NMR and other sys- 
tems [U |26] . To understand the underlying physics more 
clearly, we choose the interaction picture and make a 
transformation to a rotating frame, denoting the ground 
state and excited state of the qubit (TLS) as |0) and 
|1) (\g) and |e)), respectively. In the resonant case, i.e., 
u) = Lo q = wt, the Hamiltonian for the coupled system 
can be simplified as [TT] 
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where the basis states are \0g), |lg), |0e), and | le) (in- 
set of Fig. 1(a)). It is noticed that the Hamiltonian 
([I]) has the same form as that of the four coupled quan- 
tum pendulums |27j . The simplest way to analyze such 
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system is to find the stationary solutions without consid- 
ering dissipation, which only affects the amplitude of the 
probability. The eigenvalues of the Hamiltonian ([!]) can 
be easily obtained in the form: Ai = (O mc + f2 c )/2, A2 = 
(^mc-^c)/2, A3 = — (O mc — O e )/2, A4 = — (Q mc +fi c )/2, 
with Q mc — \/f2^ + Og. The time evolution of the 

system can be described by \^(t)) = Ylt=i c i\ipi)e~ lXit , 
where is the eigenfunction of the Hamiltonian (1) 

corresponding to A,-. Thus we can obtain the probability 
of being in state \<j>) {<j> = Og, lg, Oe, le): 



P<h = 



c*Ci(ipj\4i){^\ijji)e 



-i(A«-Aj)i 



(2) 



Note that the temporal oscillation of P^ is composed of 
Cf = 6 frequencies. However, only four frequencies are 
observable because there are two pairs of double degen- 
eracies: fii = |Ai — A4I = Cl mc + O c , VI2 = I A2 — A3 1 = 
^mc — ^c: ^3 = l^i — A3 1 = I A2 — A 4 | = Q mc , and 
= |Ai — A2I = I A 3 — A4I = fi c . Assuming the sys- 
tem is initially prepared in 1^(0)) = \0g), it is easy to 
obtain P^: 

Pog = i wr[ ! £$w + ^ cos fiit + ^ cos Q 2 t + cos fi 3 i 

+ cos n^t] 

P^s = 4a ^2 [1 — I cos cos U 2 t — cos Q^t + cos Hit] 

Poe = 4a N 2 [1 + 1 cos Qit + I COS fl 2 t — COS Q^t — COS f^t] 

p _ 4a 2 b 2 r q 4 + 6 4 a 2 n,f 
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where a = fl m /2, b — (O mc + f) c )/2 and iV = 
[f2^ + (Q mc -I- f2 c ) 2 ] /2 being the normalization factor. 
Apparently, the populations will oscillate in time with 
four frequencies. These analytical results form the foun- 
dation for a thorough understanding of the coupled qubit- 
TLS system. Moreover, these general results are also 
valid for a wide variety of four-level quantum bipartite 
systems with resonant ac drive. Based on the results 
many seemingly counterintuitive phenomena in this type 
of four-level systems become straightforward to under- 
stand both qualitatively and quantitatively. 

In the experiments involving the TLSs, the quantity 
one usually measures is the probability of finding the 
qubit in state |1), i.e., Pi = P lg + P\ e : 



Pi = 
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It is interesting to notice that although there are four 
frequency components in each P^ (Fig. 2(a) and Fig. 
2(b)), only two of them appear in Pi (Fig. 2(d)). More- 
over, fl m and £l c define three regimes for the bipartite 
system, which show significantly different behaviors: (i) 
The strong field regime, Q m /Ct c ^> 1; (ii) The interme- 
diate field regime [H [25], f2 m /f2 c ps 1 ; (iii) The weak 
field regime [21], £l m /Vl c <C 1. It should be pointed out 
that previous works focused only on one of the regimes 
thereby they could not capture the complete picture. 
Below, based on Eq. (4), we discuss the behaviors of 



the resonantly driven qubit-TLS system in the three dif- 
ferent regimes one by one. The predicted phenomena 
are demonstrated experimentally by measuring the spec- 
troscopy and the coherent oscillations in a superconduct- 
ing phase qubit coupled to the TLSs. The detailed ex- 
perimental setup and procedures have been described 
elsewhere |16| . Shown in Fig. 1 are examples of the 
spectroscopy and the coherent oscillations. The splitting 
caused by the qubit-TLS interaction is clearly observed 
at 16.572 GHz giving the coupling strength Sl c /27r 26.5 
MHz. Away from the splitting Rabi oscillation induced 
by the microwave field has the usual damped sinusoidal 
form (Fig. 1(b)). When the qubit is biased at the split- 
ting, anomalous oscillations are observed, with interest- 
ing features determined by the microwave amplitude as 
described below. 

(i) Strong field regime: Rabi beating 

In the strong field limit, i.e., fi m /f2 c > 1, 5 ra o = 
Q m /2, Pi has a simple form: 

Pi = -[2 — cos(fi TO + Q c )t — cos(fi m — il c )t}. (5) 

It is known in acoustics that beating happens as an inter- 
ference between two waves of slightly different frequen- 
cies. Here the two frequencies in Eq. ^ are close to each 
other, i.e, il m + il c ~ tt m — f2 c « f2 m , which satisfies the 
condition of beating very well. To be more clear, we fur- 
ther write Eq. ^ as Pi = |(1 — cosfi TO icosfi c i). This 
is demonstrated experimentally in Fig. 1(c), in which Pi 
appears to oscillate at f2 m with the amplitude modulated 
by a much lower frequency tt c . 

To quantitatively characterize the Rabi beating de- 
scribed above, we define Qb = Q m /2Q C as the frequency 
contrast of the beating, which is the number of Rabi oscil- 
lation periods between the two nodes of the slow varying 
envelope. Shown in Fig. 2(c) are the measured oscil- 
lations of Pi at various microwave powers. As the mi- 
crowave power is increased, Qb increases from about 2 
to 7. Therefore the beating becomes increasingly clear. 
The largest Qb is determined by Q m , Q c , and the de- 
coherence time. The exponential decay of the beating 
envelop is due to the energy relaxation. In the previous 
experiments [H [22] |2S] , anomalous Rabi oscillations due 
to the coupling between two qubits and qubit-TLS were 
reported. However, since the condition f2 m 3> Q c is not 
fulfilled, Qb is generally less than 2 and thus no clear 
pattern of Rabi beating was observed, although several 
theoretical works 23 25, 29 have predicted its existence 
in the superconducting qubits. Nevertheless, one can al- 
ways apply the Fourier transform (FT) to obtain the two 
frequency components, as shown in Fig. 2(d). With 
the help of FT, low-Qf, Rabi beatings can be revealed. 
Therefore, we argue that our experiment is the first to 
clearly demonstrate Rabi beating in a superconducting 
phase qubit. Quantum beating is usually found among 
the three-level atomic systems and has been applied to re- 
solve the detailed internal structure of matter [3U] . Rabi 
beating in the coupled qubit-TLS systems thus can be 
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utilized as a powerful tool in investigating the origin and 
properties of the TLS. 

In Fig. 2(d), we find that the difference between 
the two Rabi frequencies obtained from FT is exactly 
2f2 c , which is independent of the microwave power as 
expected from Eq. (4). In addition, it is interesting 
to notice that in the strong field limit the population 
of finding the TLS in the state \e) has a simple form 
P e = P e + Pie = (1 — cos O c t) /2, which can be viewed as 
Rabi oscillation between the subspaces {|0g), |lg)} and 
{|0e),|le)} [IT]- The oscillation frequency (f2 c ) between 
these two subspaces is half of that (2il c ) between \lg) 
and |0e) since the probability of finding the system in 
\lg) (|0e)) in each subspace is exactly 1/2. 

(ii) Intermediate field regime: anomalous Rabi oscilla- 
tion 

In this regime, f2 m ~ Q c , the two frequencies in Pi are 
well separated. The frequency contrast of the beating 
Qb is close to unity and no clear pattern of Rabi beating 
can be observed, as shown in Fig. 2(c). In particular, in 
the region of weak ac driving interplay between the two 
frequencies results in anomalous Rabi oscillations which 
have been observed previously in various superconduct- 
ing qubits [H[M1[5S]. We can extract O c and fi m by using 
FT of the anomalous Rabi oscillations. In this regime, 
the qubit is most strongly affected by the TLSs. Thus, 
much care has to be taken when performing quantum 
information processing if the Josephson junction is the 
intended qubit only. Since b is always greater than a, the 
weight of the fl\ component is smaller than that of f^2 
from Eq. (4) (indicated by the color in Fig. 2). With fl m 
further decreasing to the weak field regime, the frequency 
fii disappears, and Pi oscillates with a single frequency 
O2 , which will be discussed in detail below. 

(iii) Weak field regime: extraordinary two-photon tran- 
sitions 



In this case, Vt m -C fl c , Pq c 



P19 = 



8a 2 b 2 



cosil 3 i) <C 1, the populations in the states |1<7) and |0e) 
are very small. The system mainly evolves in the sub- 
space spanned by \0g) and |le). This can be seen more 
clearly by checking the eigen-wavefunction of the Hamil- 
tonian ([I]) directly. In the weak field limit, a <C 6, we 
have 



(6) 



Obviously, starting from the initial state \0g) the sys- 
tem will evolve in the subspace spanned by \0g) and |le) 
in the form of |*(t)) = (|-02> >e~ lX2t ~ l^e"^ 3 *)/^ 





[{e-' lX2t + e- lX3t )\0g) + (e" 
easy to obtain 



-i\ 3 t 



)|le)] /2 and it is 



Pi * Pie = \{l ~ COstt 2 i) « 1(1 - COS ^t). (7) 



The oscillatory behavior is quite normal (i.e., sinusoidal) 
except the square dependence of the Rabi frequency on 
the microwave amplitude, which is the signature of two- 
photon transitions. Therefore, when the microwave field 
is weak, we expect two-photon Rabi oscillations to occur, 
as reported in [21]. Furthermore, Eq. (7) predicts that 
the stronger the qubit-TLS coupling strength is (larger 
ftc), the smaller the oscillation frequency becomes (Fig. 
2(e)). Notice that both of these results are counterin- 
tuitive because one usually would expect that a stronger 
qubit-TLS coupling would lead to a faster oscillation and 
that two-photon transitions would be significant in strong 
ac fields. These intriguing phenomena can be under- 
stood with the help of the energy level structure shown 
in the lower inset of Fig. 3. When a microwave with 2uj q 
matching the energy difference between \0g) and |le) is 
applied, population can be transferred from \0g) to |le) 
with the help of the intermediate states ^(|lp) ± |0e)), 
although there is no direct coupling between \0g) and 
I le) . However, both of the two intermediate states are 
off-resonant with the microwave field. Thus the larger 
O c , the greater the detuning, and the smaller the fre- 
quency of the two-photon Rabi oscillations [31]. Since 
the decoherence time of our qubit is relatively short, we 
use the spectroscopy data to demonstrate our prediction. 
The spectroscopy data in Fig. 3 were obtained using long 
microwave pulses. The stationary population generated 
by the microwave induced transitions is measured. No- 
tice that in the range of frequencies measured two split- 
tings resulting from qubit-TLS coupling can be clearly 
observed [IB], with 2fi c /2?r about 20 MHz at 16.590 GHz 
and 64 MHz at 16.510 GHz, respectively. Inside each 
splitting there is a stripe in the avoided crossing. These 
stripes are resonant peaks generated by the two-photon 
transitions because the positions of the peaks (in fre- 
quency) match exactly to one half of the energy difference 
between \0g) and | le) in the entire flux bias range. The 
stripe inside the 20 MHz splitting has a higher intensity 
than that inside the 64 MHz splitting, which agrees with 
our prediction as discussed above. We calculated the sta- 
tionary population of the driven qubit-TLS system using 
the Markovian master equation. Shown in the upper in- 
set of Fig. 3 is the dependence of the spectrum intensity 
on fi c . The peak induced by the two-photon transitions 
(the middle one) disappears as f2 c increases. On the other 
hand, the two-photon transition gradually vanishes while 
the single-photon transition becomes dominant as the mi- 
crowave power increases. Further increasing microwave 
power leads to the anomalous Rabi oscillations and Rabi 
beating. 

In summary, we present a theoretical model to describe 
the quantum dynamics of a resonantly driven supercon- 
ducting qubit coupled to a TLS. The analytical result 
gives a clear physical picture of the system's dynamical 
behavior and predicts that in a four-level bipartite quan- 
tum system depending on the relative strength of the 
resonant ac driving and the interparticle coupling, high- 
contrast Rabi beating, anomalous Rabi oscillations, and 
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extraordinary two-photon transitions can occur. All of 
these phenomena have been unambiguously observed in 
our experiment using a superconducting phase qubit cou- 
pled to the TLSs and the data agree remarkably well with 
the quantitative prediction of the model. We emphasize 
that our model not only provides a unified theoretical de- 
scription of and physical insights into various experimen- 
tal observations in the superconducting qubits reported 
in the literature, but can also be applied to understand 
the dynamics of many other resonantly driven four-level 
quantum bipartite systems. The model thus forms a solid 
theoretical foundation and provides clear physical intu- 
ition to the design and analysis of coupled bipartite qubit 



systems for the quantum information processing. 



Acknowledgments 

This work is partially supported by NCET, NSFC 
(10704034, 10725415), the State Key Program for Basic 
Research of China (2006CB921801) and NSF Grant No. 
DMR-0325551. We acknowledge Northrop Grumman ES 
in Baltimore MD for technical and foundry support and 
thank R. Lewis, A. Pesetski, E. Folk, and J. Talvacchio 
for technical assistance. 



[1] Y. Makhlin, G. Schon, and A. Shnirman, Rev. Mod. 

Phys. 73, 357 (2001). 
[2] J. You and F. Nori, Physics Today 58, 42 (2005). 
[3] J. Clarke and F. K. Wilhelm, Nature 453, 1031 (2008). 
[4] R. W. Simmonds, K. M. Lang, D. A. Hite, S. Nam, D. P. 

Pappas, and J. M. Martinis, Phys. Rev. Lett. 93, 077003 

(2004). 

[5] J. M. Martinis, K. B. Cooper, R. McDermott, M. Stef- 

fen, M. Ansmann, K. D. Osborn, K. Cicak, S. Oh, D. 

P. Pappas, R. W. Simmonds, and C. C. Yu, Phys. Rev. 

Lett. 95, 210503 (2005). 
[6] Z. Kim, V. Zaretskey, Y. Yoon, J. F. Schneiderman, M. 

D. Shaw, P. M. Echternach, F. C. Wellstood, and B. S. 

Palmer, Phys. Rev. B 78, 144506 (2008). 
[7] I. Martin, L. Bulaevskii, and A. Shnirman, Phys. Rev. 

Lett. 95, 127002 (2005). 
[8] L. Faoro and L. B. Ioffe, Phys. Rev. Lett. 96, 047001 

(2006). 

[9] M. Constantin, C. C. Yu, and J. M. Martinis, Phys. Rev. 

B 79, 094520 (2009). 
[10] C. Miiller, A. Shnirman, and Y. Makhlin, Phys. Rev. B 

80, 134517 (2009). 
[11] J. M. Martinis, Quantum Inf. Process. 8, 81 (2009). 
[12] R. W. Simmonds, M. S. Allman, F. Altomare, K. D. 

O. K. Cicak, J. A. Park, M. Sillanp, A. Sirois, J. A. 

Strong, and J. D. Whittaker, Quantum Inf. Process. 8, 

117 (2009). 

[13] A. M. Zagoskin, S. Ashhab, J. R. Johansson, and F. Nori, 

Phys. Rev. Lett. 97, 077001 (2006). 
[14] M. Neeley, M. Ansmann, R. C. Bialczak, M. Hofheinz, E. 

L. N. Katz, A. OConnell, H. Wang, A. N. Cleland, and 

J. M. Martinis, Nature Physics 4, 523 (2008). 
[15] Y. Yu, S.-L. Zhu, G. Sun, X. Wen, N. Dong, J. Chen, P. 

Wu, and S. Han, Phys. Rev. Lett. 101, 157001 (2008). 
[16] G. Sun, X. Wen, B. Mao, J. Chen, Y. Yu, P. Wu, and S. 

Han, Nature Communications 1:51 (2010). 
[17] X. Wen, S.-L. Zhu, and Y. Yu, Phys. Rev. B 80, 094507 

(2009). 

[18] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
Cambridge, 2000). 

[19] F. Deppe, M. Mariantoni, E. P. Menzel, A. Marx, S. 
Saito, K. Kakuyanagi, H. Tanaka, T. Meno, K. Semba, 
H. Takayanagi, E. Solano and R. Gross, Nature Physics 
4, 686 (2008). 



[20] D. Hanneke, J. P. Home, J. D. Jost, J. M. Amini, D. 
Leibfried, and D. J. Wineland, Nature Physics 6, 13 
(2010). 

[21] A. Lupascu, P. Bertet, E. F. C. Driessen, C. J. P. M. Har- 

mans, and J. E. Mooij, Phys. Rev. B 80, 172506 (2009). 
[22] J. Lisenfeld, C. Miiller, J. H. Cole, P. Bushev, A. 

Lukashenko, A. Shnirman, and A. V. Ustinov, Phys. Rev. 

B 81, 100511 (2010). 
[23] Y. M. Galperin, D. V. Shantsev, J. Bergli, and B. L. 

Altshuler, Europhys. Lett. 71, 21 (2005). 
[24] L.-C. Ku and C. C. Yu, Phys. Rev. B 72, 024526 (2005). 
[25] S. Ashhab, J. R. Johansson, and F. Nori, New J. Phys. 

8, 103 (2006). 

[26] L. M. K. Vandersypen and I. L. Chuang, Rev. Mod. Phys. 

76, 1037 (2005). 
[27] R. P. Feynman, The Feynman Lectures on Physics III 

(Addison- Wesley, 1989). 
[28] Y. A. Pashkin, T. Yamamoto, O. Astafiev, Y. Nakamura, 

D. V. Averin, and J. S. Tsais, Nature 421, 823 (2003). 
[29] J. Claudon, A. Zazunov, F. W. J. Hekking, and O. Buis- 

son, Phys. Rev. B 78, 184503 (2008). 
[30] D. Budker, W. Gawlik, D. F. Kimball, S. M. Rochester, 

V. V. Yashchuk, and A. Weis, Rev. Mod. Phys. 74, 1153 

(2002). 

[31] A. F. Linskens, I. Holleman, N. Dam, and J. Reuss, Phys. 
Rev. A 54, 4854 (1996). 



5 




FIG. 1: Spectroscopy and coherent oscillations, (a) Spec- 
troscopy of the qubit versus the flux bias with a splitting at 
/ = 16.572 GHz due to the coupling of the qubit-TLS. (b) 
Usual Rabi oscillation with the damping time Tr as 81.5 ns 
at / = 16.728 GHz (arrow in the left) where the effect of 
the TLS is negligible, (c) Coherent oscillation at the avoided 
crossing (arrow in the right) shows quantum beating due to 
the interference of Rabi oscillations in the coupled system. In 
(b) and (c), the red dots are the experimental results and the 
solid lines are the theoretical results. 
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FIG. 2: Frequencies Qr in P^. (a) and (b), Four frequencies 
with different weight (indicated by the color) in P\ a and P\ e 
versus Q m , respectively, (c) Rabi oscillations with the mi- 
crowave power, at the top of the fridge, increasing from -13 
dBm to -1 dBm with a step of 1 dBm from bottom to top. 
Curves are shifted vertically for clarity. Quantum beating be- 
comes more clear as the amplitude of microwave increases, 
(d) Frequencies (dots), obtained by the Fourier transforma- 
tions of the corresponding Rabi oscillations in (c), versus the 
microwave amplitude. The color lines are the two frequencies 
in Pi obtained from the theoretical analysis, (e) Frequencies 
in Pi induced by the two-photon transitions versus Q m with 
three different Q c /2tt: 15.0 MHz, 26.5 MHz, and 50.0 MHz. 



G 



Ol 02 03 

i : : : 

■0.5 




8.18 8.20 8.22 8.24 8.26 8.28 

Flux Bias (arb. units) 

FIG. 3: Two-photon transitions. Large f2 c hinders two- 
photon transitions, leading to a lower resonant peak in the 
measured spectrum as marked with arrows. The lower in- 
set shows the energy levels in the qubit-TLS coupled system. 
The upper inset shows the spectrum density versus fi c . The 
outside peaks are due to the stationary population from |0<;} 
to ^Tjdlfl) ± |0e}). The middle peak is due to the two-photon 
transitions from \0g) to |le). 



